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EFFECT OF COMPRESSIBILITY AT HIGH SUBSONIC VELOCITIES ON THE LIFTING FORCE 

ACTING ON AN ELLIPTIC CYLINDER 


By Carl Kaplan 


SUMMARY 

An extended form oj tke Ackeret iteration method, applicable 
to arbitrary profiles , is utilized to calculate the compressible flow 
at high subsonic velocities past an elliptic cylinder. Tke angle 
i rf attack to tke direction of the undisturbed stream is small and 
ike circulation is fixed by the Kutta condition at the trailing end 
oj the major axis. The expression jor tke lifting force on the 
elliptic cylinder is dented and shows a first-step improvement 
of tke Prandtl-Glauert rule. It is further shown that the 
expression for tke lifting force, although derived specifically for 
an elliptic cylinder, may be extended to arbitrary symmetrical 
profiles. 

INTRODUCTION 

Two methods of approximation, complementary in nature, 
have been used for the solution of problems of flow past 
bodies in the subsonic-velocity range. They are the Janzen- 
Rayleigh method (reference I), in which the velocity po- 
tential or the stream function is developed in a power series 
of the stream Mach number, and the Ackeret iteration 
method (reference 2), in which the velocity potential or the 
stream function is developed in a power series of a geometrical 
parameter characteristic of the body. The complementary 
nature of these two methods lies in the fact that the Janzen- 
Rayleigh procedure yields accurate results in the case of 
thick bodies, for which the critical stream Mach numbers are 
low; whereas the Ackeret iteration process yields accurate 
results in the case of slender shapes, for which the critical 
stream Mach numbers are in the neighborhood of unity. 

In the Ackeret iteration process the assumption is made 
that, if X is a parameter that represents the departure of the 
profile shape from a straight-line segment at zero angle of 
attack (that is, X may be the thickness, the camber, or the 
angle of attack), the stream function j/, say, may be developed 
in a power series of X in which the coefficients are functions 
of the plane coordinates x and y and of the stream Mach 
number d/ r ; that is, * 

4>= — Uy+ft(x, y; Aft'fi+ffx, y; M t )X*- f- ... (1) 

where U is the velocity of the undisturbed flow. By insert- 
ing this expression into the fundamental nonlinear differ- 
ential equation for f and by equating the coefficients of the 
various powers of X to zero, a system of linear differential 
equations for the functions f n {x,y;M,) is obtained. The 
integration of this system of differential equations can be 
performed for as many steps as desired, the first step 


ft (x,y;M,) being equivalent to the Prandtl-Glauert approxi- 
mation. There is, however, a fundamental difficulty with 
the form of the development, equation (1), which does not 
appear in the Janzen-Rayleigh method. In the Janzen- 
Rayleigh method the expansion of the stream function in 
powers of the stream Mach number, namely, 

i'lAIf -f- . . . (2) 

can always be obtained, whereas the possibility of the ex- 
pansion, equation (1), cannot be guaranteed beforehand. 
This difficulty in the Ackeret process was avoided in refer- 
ences 3 and 4 by choosing as solid boundaries profiles having 
no stagnation points. In such cases, a development of the 
form of equation (1) is always possible. When, however, 
shapes are chosen that possess stagnation points, terms of 
the type X* log X ultimately appear on the right-hand side 
of equation (1) and the explicit development in powers of X 
is strictly no longer possible. This difficulty may be avoided 
by assuming a somewhat more general form for the develop- 
ment of the stream function; namely, 

t=-Uy+Mx, y)+M*, V)+ • ■ • (3) 

where the shape parameter X is contained implicitly in the 
various functions f n - In equation (3) fa corresponds to the 
Prandtl-Glauert approximation and 4/ is made to satisfy the 
exact boundary conditions at the solid and at infinity. For 
the purpose of defining and controlling the iteration procedure, 
the function ^ B+ i is regarded as small compared with the 
preceding function fa and the derivatives have a similar 
relationship. It can be stated now that the aforementioned 
difficulty is more apparent than real, for in carrying out the 
iterative steps according to equation (1) X behaves like a 
dummy parameter, which serves only to regulate the itera- 
tion process in exactly the same manner as equation (3). 

For slender bodies the first few steps of the Ackeret itera- 
tion method may be expected to yield a good result with the 
exception of a small region in the neighborhood of a stagna- 
tion point Hantzsche, in reference 5, treated the case of an 
elliptic cylinder at zero incidence in a uniform stream accord- 
ing to Ackeret’s process as represented by equation (1), with 
the exact boundary conditions at the solid and at infinity 
being satisfied. These calculations were carried through 
the X*-t-erms and included a term X 4 log X. A comparison, 
made in reference 5, of this result with that obtained by the 
present author according to the Janzen-Rayleigh method 
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showed agreement in the terms common to the two develop- 
ments. A similar comparison made for the “bump” in 
reference 3 also showed complete agreement in the terms 
common to the two methods. These comparisons illustrate 
the fact that solutions obtained by the Janzen-Rayleigh 
and the Ackeret methods are simply different representations 
of a single unique solution. In particular, it is evident that 
the Janzen-Rayleigh method, which is valid at stagnation 
points, includes parts of the Ackeret development and, 
conversely, the Ackeret development includes parts of the 
Janzen-Rayleigh result. Although the first few terms of 
equation (3) do not yield very accurate results at stagnation 
points, these terms nevertheless represent correctly, to some 
extent, the effect of compressibility at such points. The 
accuracy of the calculations obviously depends on the num- 
ber of terms derived, each new term reducing the extent 
of the region of inaccuracy in the neighborhood of a stagna- 
tion point. 

The question of the convergence of the sequence of func- 
tions, equation (3), is a difficult one and should be thoroughly 
investigated. Schmieden and Kawalki (reference 6) state 
that both the Janzen-Rayleigh and the Ackeret develop- 
ments diverge when the local velocity of sound is first ex- 
ceeded in the region of flow, that is, at the critical value of 
the stream Mach number. This statement is contradicted, 
however, by results obtained by means of hodograph or 
velocity variables. In general, as shown in reference 7, the 
limit of steady potential flow of a compressible fluid is de- 
termined by the vanishing of the Jacobian of the transforma- 
tion from physical x, y-variables to hodograph 0, g-variables. 
The vanishing of .this functional determinant does not neces- 
sarily occur when the local Mach number first readies unity 
and, consequently, there exist continuous solutions of the 
general differential equation governing the flow of a com- 
pressible fluid for which a part of the region of flow is super- 
sonic. It is reasonable to conclude, therefore, that the series 
solutions given by equations (2) and (3) diverge at the .value 
of the stream Mach number that marks the limit of poten- 
tial flow rather than at the value for which the local velocity 
of sound is first exceeded in the flow. 

The Ackeret process in the form of equation (1), but for 
the velocity potential, was used to calculate the flow past a 
bump and past a circular arc profile in references 3 and 4, 
respectively. This calculation was possible because the 
profiles considered did not possess stagnation points. As a 
consequence, the problem could be treated completely in the 
physical plane and, moreover, the boundary conditions on 
the vdocity potential were tractable. In general, however, 
when shapes with stagnation points are treated, the Ackeret 
process in the form of equation (3) instead of in the form of 
equation (1) should be utilized. In such cases the problem 
is most conveniently treated in a new plane related to the 
plane of flow by an affine transformation. In this plane, 
however, the boundary conditions for the velodty potential 
beeome very complicated. On the other hand, the boundary 
condition for the stream function, ^=0 at the solid, is 
invariant for affine transformations and it is— therefore 
suggested that the stream function be utilized. The trans- 


ference. of the problem to an affinely connected plane un- 
fortunately introduces a distortion in the actual profile that 
varies with the stream Mach number. In the case of an 
elliptic profile, however, this distortion does not matter, for 
the affine distortion of an ellipse leads to another ellipse. 
This property of the ellipse makes it a preferred profile for 
the Ackeret iteration process in the form of equation (3) and 
is also the reason it is chosen as the solid boundary for the 
problem treated in the present paper. Specifically, the prob- 
lem treated herein is the determination of the flow r of a com- 
pressible fluid past an elliptic cylinder set at a small angle of 
incidence in a uniform stream, with circulation determined 
by the Kutt-a condition at the trailing end of the major axis. 
The main purpose of this calculation is to obtain some infor- 
mation with regard to the effect of compressibility at high 
subsonic stream Mach numbers on the lifting force acting 
on an elliptic cylinder. A calculation is now in progress in 
the Physical Research Division of the Langley Memorial 
Aeronautical Laboratory to find the effect of compressibility 
at high subsonic stream Mach numbers on the moment and 
on the location of the center of pressure of an elliptic cylinder. 

CALCULATION OF THE FIRST AND SECOND APPROXIMATIONS 
The equation of continuity can be written as 

and the condition for irrotational flou r as 


bp bu 


( 5 ) 


where. 

X, Y rectangular Cartesian coordinates in physical flow plane 
u, v components of velocity along A r -axis and I r -axis, 
respectively 

p variable density of fluid 

Pi constant density of undisturbed fluid at infinity - 


Equation (5) defines a velocity potential <j>, where 


_b£. 



( 6 ) 


and equation (4) defines a mass flow or stream function 
where. . 


_d<£ _Pid\p \ 

~dX~ P W 

— _Pi 
~bF p 




( 7 ) 


From equations (4) and (5), with the use of equations (6) and 
(7), respectively, the following alternate forms of the basic 
differential equation are obtained: 


and 


5 /p b<£\ j_ 5 ( p b<£\ 
b />! b^ \ , b /pi bA_ 

bzV7W + 5FV7aiy~ 0 


( 8 ) 

( 9 ) 
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For reasons stated in the “Introduction,” equation (9) for 
the stream function will be treated in the present paper. 
The difficulty of obtaining a solution lies in the fact that the 
density of the fluid p is related to u and v because of the 
dependence of p on the local pressure. If it is postulated 
that the fluid is nonviscous and if the fluid is assumed to 
flow from infinity with a constant velocity U, the total energy 
of the fluid will then have the same value at every point in 
the region of flow. If p and q denote the pressure and the 
velocity of the fluid, respectively, Bernoulli’s equation 
becomes . 

( 10 ) 

where the lower limit of integration refers to the starting 
conditions at infinity. Moreover, p and p are connected 
by some known adiabatic equation of state such as 


Constant (11) 

p v Pi t 

where, in the case of a perfect gas, 7 denotes the ratio of 
specific heats at constant pressure and at constant volume. 
By means of the definition of the velocity of sound, 


the following relations are obtained from equation (10) : 

c*=ci s [i-^r M * (&- 0] (*») 

P=Pi[ 1 -^M 1 , (^- 1 )]- (14) 

P=Pi [l • W (jjs- l)J^ (15) 


where M t = — is the Mach number of the undisturbed stream 
at infinity and q*=u , +v l . 

By means of equations (7), equation (14) can be rewritten 
in the form 


H 1 M ' &' 0+7 - 0 


where p x and p Y denote dp/dX and tipjdY, respectively. If, 

for the moment, pi/p and - x ^72 ~ 1 are considered to be 

dependent and independent variables, respectively, a 
Maclaurin expansion in the neighborhood of the undisturbed 

t 2 1 1 2 

stream, where * x —1=0 and ^ = 1 , yields 


a=i+±<7~i>(*^'-i) 

+|(m 2 -i)M(y+4)+(7+D (^-1)3 • • ■ 

(17) 

S42951 — 49 6 


where 2 and ju 2 — The form of this 

development has been chosen to be consistent with the 
Ackeret iteration process, which is essentially an iteration 
around an undisturbed stream of zero incidence. Corre- 
sponding to equation (3), 

P= — UY-\-pi-\-p a -\- • . . (18) 

where p n+1 is regarded as small compared with p K and the 
derivatives have a similar relationship. When this ex- 
pression for p is substituted into equation *(17) and 1 p ir is 
noted to be of the same order as 1 p^ or Piy 1 , then inclusive 
of the second power in the derivatives, 

7=1-^ i) 7 _1) \yt — 2 lP~) 

+|(M 2 -l) 2 [(7+4) + (7+D(M a -l)]%- 2 + . . . (19) 

When this expression for pi/p is substituted into the basic 
differential equation (9) and terms of equal order of magnitude 
in the derivatives of p n are equated, the following differential 
equations for ip\ and p 2 are obtained: 

TplXX-\- l^PlYY—O ( 20 ) 

and 

pixx^ MV 2 rr= (p 2 — 1) Pixx~\r 2 'Pixy 

+M 2 [3+(7+l)(M 2 -D]7f^irr} (21) 

These differential equations are most easily solved by first 
applying the affine transformation 

2= AT 

y=S Y 

Equation (20) then becomes a Laplace equation 

Plxx+Plvi^O (23) 

and equation (21) becomes a Poisson equation 

Pixx~\~p2w = {f^ 1) J^JJPlxX+2 ~^Jj'l'lX¥ 

+ [2+{y+W-D]fyPi„} (24) 

For purposes of calculation a new stream function p* may 
be conveniently introduced, where 

P=ixUp* (25) 

From equations (22) it can be seen that 

Px=Px=nUp* s 

and 

■ Pt= Up* v 

so that the undisturbed flows, at zero angles of incidence, 
in the physical and affinely distorted planes are identical; 
that is, (Py)i=(UP*v)i- With the introduction of the 
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stream function $*, equations (23) and (24) become 

= 0 (26) 

and 

— l){'P*it'P*ixx- l r2'l , *ix4'*ixi 

+ [3 + (Y+l)(/x i -l)]^ 1 ^*i„} 

or, with the use of equation (26), 

^**„+^** w =(M*-l){-[2+(Y+l)0* , -l)]f*n^* 1 .. 

. +2** lx t* lxv } - (27) 


Equations (26) and (27) show that, in order to calculate 
the various approximations, the incompressible flow past 
the distorted profile in the zy-plane must be known. Thus, 
if a profile is given in the physical XT-plane, it will be neces- 
sary to find the conformal transformation to a circle of the 
distorted profile in the zy~plane for each value of the stream 
Mach number, since the affine distortion of a profile depends 
on the stream Mach number, as shown by equation (22). 
In general, then, the problem to be solved is the flow past 
an arbitrary profile in tho affinely distorted ry-plane. This 
procedure in general involves the laborious calculations of 
the coefficients of conformal transformations to a circle for a 
number of values of the stream Mach number in such a 
way that the distorted profiles in the ry-plane correspond 
to the given profile in the physical XT-plane. In the case 
of an ellipse, however, the distorted profile is again an ellipse, 
and it is therefore a simple matter to transfer the results 
obtained in the zy-plane to the XT-plane of the original 
elliptic profile. In the present paper the elliptic profile is 
so orientod that its major axis lies along the X-axis. The 
relation between the profiles in the two planes is then given 
by 


a=a' 



V- 


u 2 — 1 

b n 

b 

R=R'+^b' 

2/x 


(28) 


(Wi'+Wii) 

&zif = 2 Wjjj) 

(Witz WiJi) 

where a bar indicates conjugate-complex quantities. There- 
fore, 

2&^i*,=I. P. (w ui w u +wi„wrz) = (w„+ W]7) I. P. w lt , 

and 

fyixfox^l. P. {WutWu — Wi„u'u) =—i(w u —Wu) R. P. u'u, 

Equation (27) can then be written as 

G« 2 -l) I.P. |j <Kt*. , ).-KH-4)(«* 1 W,;) I ] (29) 
where <r=( 7 -[-l)(/t 2 — 1). 

Consider now a circle of radius R, with its center at the 
origin of the z'-plane, into which the distorted profile in 
the z-plane is mapped by means of a conformal transforma- 
tion. Any point on the circle can be expressed in the form 

z'=x' +iy'=R (cos £— i sin £)=2?e -< f 

so that the point describes the circle in the clockwise sense 
and leaves the region outside the circle on the left. If now 
t=Z+iri, the transformation 

z'—Rer't (30) 

yields a circle of radius R when t?= 0 and the infinite region 
of tho z'-plane when jj= 4- <». Equation (30) can be looked 
upon as the transformation from Cartesian coordinates x‘ , y' 
to polar coordinates Re*,—t. The conformal transformation 
of an arbitrary profile in the z-plane into a circle of radius R 
in the z'-plane, with center at the origin, can therefore be 
written 

*=/(«' 0 (31) 


tan a=- tan a' (or, to the first order, a = — «') 


where 


R — 2 (tt+fi) 


a semimajor axis of ellipse in affinely distorted plane 
b semiminor axis of ellipse in affinely distorted plane 
c semifocal distance of ellipse in affinely distorted plane 
a angle of incidence of ellipse in affinely distorted plane 

and the prime indicates corresponding values in the physical 
plane. 

The solution of equation (26) is in general the imaginary 
part of an analytic function vjj (z) where z=x\iy. (The 
asterisk has been dropped.) It is easy then to verify that 


where _rj=0 yields the parametric equations of the profile in 
the z-plane. 

When f and p are introduced as independent variables, 
equation (29) becomes 

||p-|(M i -l)l.P.[|v(u^ 2 ) f ^+(«r-|-4)(u. 1 0 r tCLf](32) 

The general solution of equation (32) or of equation (29) 
is obtained directly and is 

U*-l) I.P. *zu> u 2 -K<r+4) w lWlt +^(I)] (33) 

where F(£) is an arbitrary function to be determined accord- 
ing to the boundary conditions at tho solid surface and at 
infinity. 
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Consider now an ellipse in the affine 2 -plane with semi- 
major axis a and semiminor axis b. The undisturbed stream 
at infinity makes a small angle a with the negative direction 
of the x-axis, and the circulation is determined according to 
the Kutta condition that the downstream end of the major 
axis is a stagnation point. In accordance with the Ackeret 
process, 


1 +^ 2 + • • • (34) 

so that the stream function is expanded around an undis- 
turbed stream at zero incidence. Since the angle of attack 
has been assumed small, powers of a higher than the first 
are neglected. The boundary conditions are then, at the 
surface of the ellipse, ij= 0 , 

fa =V 
fa =0 


and, at infinity, ij= a> , 


bfa 

dr c 


(35a) 


bfa 

by 


=o 



(35b) 



Now, the transformation 

(36) 

maps the region external to the ellipse in the 2 -plane into 
the region external to a circle of radius R=]i (a + 6 ) with 
center at the origin of the z'-plane. The complex potential 
for the flow considered is given by 

w=—(z'e ia +^-^f — 2fi? sin a log ^ (37) 

When the variable f is introduced by means of equation (30), 
equation (36) becomes 

z—c cos (£+iX) (38a) 

and equation (37) becomes 

W——2R cos (£— a)— 2 .Rf sin a (38b) 

where a=c cosh X and b=e sinh X. This expression for w 
includes the uniform undisturbed stream —z which must be 
extracted in order to obtain fa. Thus, for a small angle of 
attack a, 

w^=c cos (£-f-iX)— 2R cos £— 2 R (sin f+£)a: (39) 


and 


i^i=I. P.i0i 

= — c sinh 0?+X) sin £+2 R sinh q sin £—2 R (sinh q cos £+ij)a 


It can be easily verified that this expression for fa satisfies the boundary conditions stated in equations (35) and also 
that the downstream end of the major axis is a stagnation point. 

Similarly, the most direct way to determine fa is to consider it to be the imaginary part of a nonanalytic function 
of £ and f . Thus, from equations (33), (38), and (39), 


w t =— g (/i s - 1) 


^ <rc cos (?— iX)i 0 lJ i +(ff+ 4 )c[cos (?— iX)- 2 j cos £—2 j a (sin ?)]i 0 i«+F(£) 


In order to satisfy the boundary condition fa=0 at the surface, q=0, it is a simple matter to supply the functions 
of £ needed to make the coefficients of w u and w u 2 vanish for 13 = 0 . For example, cos (£— iX)= cos (£— iX) for tj= 0. 
Thus, 


fa=-±(S-l)I.J>. 


2 ffc [cos (J— iX)— cos (f — iX)]Wi , 2 


+ (<r+4)c cos (f— iX)— 2 ^ cos f—2 a(sin ?+£)— cos (f — iX)+2 ^ cos £+2 ^ a(sin £+£) kxJ 

— C C C C _ J 


(40) 


In order to satisfy the boundary condition at infinity and the condition that the downstream end of the major axis 
be a stagnation point, the procedure, according to equations (7), (22), and (25), is as follows: 


P f 1 \ n . cff d * 2 d£ d . _ . 

TJJ V Tz 5£ = Ji b£ (^-^) 

From equations (39) to (41) it then follows that at infinity 

I ( ^ _1) (ff+4) ZTb ia 


(41) 
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Hence, in order to satisfy the boundary conditions at infinity, 0 and r a =— ^=0, a term — fy*— l)(H-4)a sinf, 

the imaginary part of which vanishes for 17=0, must be added to the expression for w 2 . This addition to w t introduces 
a velocity at the downstream end of the major axis, given by 


P 

PiU 


\ P /t-r, y-0 


(P*— 1) (<r+4) ia 


Again, in order to render (£=«■, 17=0) a stagnation point, a term — i b (m 2 — 1) (H-4) aj; satisfying the boundary condi- 
tions at the solid and at infinity must be added to the expression for u> 2 . Finally then 
1/^= I. P. Vh 

= ~| (p 2 — 1) I. P. 77 <tc [cos (?— iX)— cos (f— 1\)] ifli t 2 +(a-+4) c [cos (f— i\) — 2 ^ cos f 
— 2 ^a(sinf+f)— cos(f— fX) + 2 j cos f+2 ^ a (sin f+f)J Wu+26 (<r+4)(sin f+f) a (42) 

The complete expression for 1 p, obtained from equations (34), (39), and (42), is then given by 
lA=— I- P- (2 R cos f-t-2 Ra (sin f+f)+|(p 2 — [cos (?—■ iX)— cos (f— i\)]Wu 

-h(<r+4)c [cos (?— ?X)— 2 ^ cos f— 2 ^ a (sin f+f )— cos (£— iX)+2 ~ cos f+2 j a (sin f+f)JtPi, 
+26(o-f-4)of (sin f +f)J^ (43) 

where w l is obtained from equation (39) and 


CALCULATION OF THE VELOCITY COMPONENTS u AND v 


n ■ ■» « • 

The components u and v of the velocity of the compressible fluid past the actual profile in the physical A'F-plane can be 
put into the following convenient form by means of equations (7) , (22) , and (25) : 



ctz 5 f df 


(w—w) 


(44) 


where tv is the expression on the right-hand side of equation (43) with I. P. omitted. It follows that 

^ ( u ~^ ®) = i!T(f+iXT( 2 7 [_sin f ( 1 +cos r)]+ i ( ^ 8- 1} (ff+ 4) C 1 + cos f ) ^ a 

+i(/i 2 — 1) 0- j| [w u 2 sin (f— iX)+w 1 3 2 sinG'-l-iX)3— [cos (T— iX) — cos(f— iX)]w u tP,zr 

+|fy 2 — l)(v+4){[sin(f— iX)— 2^sinf+2^ct(l + cosf)J«^+[sin (f +iX)— 2ysinf +2^ a (1-f cos 

“t” [cos (?— rX) — 2j cos?— 2^a(sinf+?) — cos fy— fX)-i-2^cosf+2 ; ^a (sin f +f)Jw uf |^ (45) 

where equations' (28) provide the correspondence between the distorted ellipse in the r-plane and the actual ellipse in the 
Z-pIane. 

Equation (45) will now be utilized to calculate the lift on an elliptic cylinder in compressible flow. For this calculation 
a control contour, which eventually is taken as infinitely large, will be applied in the usual manner. The elimination 
of the variable f in equation (45) and in the equations needed for the calculation of the lift is in general impossible and in 
the present example undesirable. Since the regions at infinity correspond in the z-. and z'-planes, it is convenient to choose 
a large circle in the s'-plane as the control contour and to effect all the calculations with $ and j ; as independent vari- 
ables. The advantage gained by this procedure is a great reduction in labor, in that on a circle rj is constant and hence 
only functions of the single variable z ' (=e" <£ ) appear. The first step is to obtain developments for pu/p x V and pvjpiU in 
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the neighborhood of infinity and then to form the combination (w — iv). This calculation, according to equation (45), 

yields the following result: • 

^(«-it,)=-( 1+.V«>+‘| [l+j („'-l)(»+4) 5 Aj][(l-rt 2 '-(l+,) I] 


— la 


d~\~ b 




(46) 


where z'=e~ Since 'Py and fa, it follows from equations (17) and (46) that 

: [ x +i V-VM sis] 

{^+1 (sr?) (/t!_ 1)(ff+4) ] [(^+ 7 ^)+^ ^ (^-74)} 


+ 


+ 


(47) 


CALCULATION OF THE LIFT 


In a compressible flow as in an incompressible flow, the lift is given by (reference 8) 

L e = Pl UT c 

where T c is the circulation round the profile and where, by definition, 


T c — (^(-wdX+®rfy)=R. P. f ( u—iv)dZ 


(48) 

(49) 


Now, from equation (22), 


Hence 


Since z=c cos (f-f-iX), 


X=x 

Y= P y 




dZ=—^-YY: c sin (f+fX)df — c sin 

For the evaluation of the line integral for r„, the control contour is a large circle in the z'-plane. Therefore, ij is a 

dz' 

constant and dt;=df=i — hence, 


dZ ~~\ [(!+/*) ( z / e i+a 2')+(l ^(7 e 2^a)] 7 


(50) 


where z'=e -,f . The desired expression for u—iv, obtained from equations (46) and (47), is 

i(il-i») = -l-i / i a + 1 i[ 1 +|(/- 1)( v+ 4) ^][ (1 -^^_ ( l+^^+ (M *- 1) (^-7)]+ . . . (51) 

Equations (49) to (51), with only terms containing the factor ds'jz' contributing to the line integral for 1%, then yield the 
following result: 


T e =4icRU P 2 a. 1+g G* 8 — 1) (ff+4) 


(52) 
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If a, b, R, and a are replaced by a', b’ , R', and a' according 
to the correspondence equations (28), then for the actual 
elliptic cylinder in the physical flow plane Z, 

r e =4irR'Ua'ii-{-2irUa , b' £(1— /i)+£ G* 2 — 1) (<r+4)J 

Since the circulation in the incompressible case, M x = 0 or 
(U=l, is 

r ( =^R'Ua! 

the ratio r c /r< is given by 

Ql— m)+| (m j — 1) (<r+4)J 
With (a' + 6') and cr=(7+l) (n 2 —l), 

0*— i)+j(y+i) (m ! -i) 2 J (53) 

where t' is the thickness coefficient b'/a' of the elliptic 
cylinder in the physical flow plane. This equation represents 
a first-step improvement of the Prandtl-Glauert approxima- 
tion and reduces to that result when t'~* 0. 

Although equation (53) has been derived specifically for an 
elliptic cylinder, it will be shown that the result can be ex- 
tended to a slender arbitrary synynetrical profile. Hantzsche 
and Wendt (reference 9) derived a similar relation for the 
case of a symmetrical Joukowski profile with a sharp trail- 
ing edge. The result of Hantzsche and Wendt may be 
written as 


[m(m— l)+£ (t+ 1) (/— l) s (54) 


Note that the function of y contained between the brackets 
is the same in equations (53) and (54). This coincidence 


suggests a correspondence between the factors 


t' , e 
l-M' and l-(-e' 


A correspondence is obtained in t he fol lowing manner: 
is well known that by means of the mapping function 


It 


Z = Z '+ W 


the circle of radius c 2 /4, with its center at the origin of the 
Z'-plane, is mapped into the line segment extending from 

£ 

z=—c to z=c, and the circle of radius ^ (1+e), with its 
c 

center at Z'=« is mapped into a symmetrical Joukowski 
profile with sharp trailing edge in the Z'-plane. Now 


V 

1+f 


Tt- c * 

R m 

2 R 


(55) 


or, with R = | (1 + e), 


1 + e - 


1 + e 


i +tr 


2(l + e) 


1+e 


Thus, to the first power in e, and the correspond- 

ence between equations (53) and (54) is established. 

In the case of an arbitrary symmetrical profile, the Theo- 
dorsen method (see reference 10) is particularly well suited 

fJ 

to obtain an expression corresponding to rxp- An essential 

feature of the potential theory of arbitrary wing profiles de- 
veloped by Theodorsen is a rapidly convergent procedure for 
obtaining the conformal transformation of the profile to a 
circle, also the radius R of the circle. The coefficient of tlio 

T^-term of this conformal transformation, denoted by c 1 / 4, 

and the radius R of the conformal circle define an ellipse 


with 


Z—c cos (£+iX) 




Then, from equation (55), 


t' 

1 +t 


7=5 (l-e- a ) 


and, therefore, for an arbitrary symmetrical profile, the for- 
mula that corresponds to equation (53) may be written 


zS =s ?!“ M+ i < ‘ 1-e (**— i)+j (y+i) (/!*— i) j J 


(56). 


Table I shdws values of the ratio LJL ir calculated by 
means of equation (53), for various values of the thickness 
coefficient t' and the stream Mach number Mi (with 7= 1.4 
for air) and figure 1 shows the corresponding graphs with Mi 
as abscissa and LJL t as ordinate. 


Langley Memorial Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., May 16, 1946. 
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TABLE I 

RATIO OF LIFTS FOR COMPRESSIBLE AND 
INCOMPRESSIBLE FLOWS 


Mi 


** 

m’— M 

+0.6Ca*-l)‘ 

L.!Li 

t'=0.05 

/'■= 0.10 


f'-OJO 

0. 10 

1.0050 

1.0101 

a 0051 

L0O53 

1.0056 

1.0067 

1.0069 

.20 

1.0206 

1.0417 

.0221 

1.0217 

1.0226 

1.0235 

1.0243 

..30 

1.0483 

1.0980 

.0566 

1.0510 

1.0634 

1.0567 

L0577 

.40 

1.0011 

1. 1005 

.1212 

1.0969 

1.1021 

1.1069 

1.1118 

.45 

I. 1108 

1.2889 

.1728 

1.12S0 

1.1365 

L 1423 

1.1486 

.GO 

1.1647 

1.3333 

.2463 

1.1664 

1.1770 

1.1867 

1.1966 

.66 

1.1974 

1.4337 

.3492 

1.2140 

1.2291 

1.2429 

1.2666 

.60 

1.2500 

1.6625 

.5023 

1-2739 

1.2967 

1.3165 

1.3337 

.66 

1.3160 

L7310 

.7376 

L3510 

1.3830 

1.4121 

1.4388 

.70 

1.4003 

1.9608 

1.1144 

1.4634 

1.5016 

1.5456 

1.5880 

.76 

I. 5119 

2.2858 

1.7660 

1.6955 

1.6716 

1.7409 

1.8046 

.80 

1.6067 

2.7778 

3.0074 

1.8099 

1. 9401 

2.0589 

2.1679 

.86 

1.89S3 

3.0037 

5.7728 

2.1732 

2.4231 

2.6613 

2.8605 

.90 

3.2942 

5.2632 

13.8737 

2.9548 

3.5654 

4.1042 

4.6064 



M, 

Figure i— Hallo of lifts for oompressibte and Incompressible fluids as function of stream 

Mach number. 





